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ABSTRACT 


This  thesis  presents  an  analysis  of  the  performance  of  a 
local/mobile  radio  communications  system  utilizing  the  Aloha 
random  access  protocol  with  variable  length  packets.  The 
capture  phenomenon  due  to  the  near/far  effect  that  enhances 
the  performance  of  the  system  is  investigated.  A  tagged  packet 
will  capture  the  base  station  if  its  signal-to-interference 
ratio  exceeds  a  threshold  y0.  Because  of  the  near/ far  effect, 
users  near  the  base  station  typically  have  a  stronger  signal 
than  those  farther  away.  A  multiple  ring  model  is  used  to 
alleviate  this  problem.  Users  in  one  ring  employ  different 
retransmission  strategy  from  those  in  other  rings.  A  shorter 
retransmission  delay  is  allocated  to  users  in  rings  farther 
from  the  base  station  than  in  those  closer  to  the  base  station 
in  order  to  achieve  approximately  the  same  average  delay 
throughout  the  network  regardless  of  location. 
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I.  INTRODUCTION 


Aloha  protocol  is  a  random  access  protocol  in  which  users 
transmit  their  packets  at  will.  Because  of  the  random  times  at 
which  users  may  decide  to  transmit,  there  is  a  chance  that  two 
or  more  users  transmit  at  overlapping  times.  This  results  in 
a  collision.  Whether  the  overlapping  covers  one  bit  or  a  whole 
packet,  the  collided  or  overlapped  packets  will  be  discarded. 
Aloha  protocol  has  inherent  characteristics  of  instability.  As 
the  traffic  intensity  increases,  the  probability  of  collision 
increases  as  well,  leading  to  possible  instabilities  of  the 
system.  Throughput  is  limited  to  some  maximum  value  less  than 
the  channe]  capacity.  The  particular  values  depend  on  the 
access  mechanism.  For  example,  the  maximum  value  for  pure 
Aloha  is  0.5e_1  or  0.18.  Despite  its  poor  throughput  and  its 
inherent  instability,  the  Aloha  protocol  is  still  a  prime 
candidate  for  many  local  and  mobile  packet  radio  systems 
because  of  its  ease  in  implementation.  In  many  practical 
applications,  Aloha  performance  in  local  and  mobile 
communication  systems  is  affected  by  the  near/far  effect  as 
well  as  multi-path  fading.  The  near/far  effect  arises  when 
users  are  located  at  different  distances  from  the  base 
station.  The  power  levels  of  the  received  packets  may  vary 
substantially  due  to  spatial  attenuation  of  the  signals,  and 
the  packet  with  highest  power  level  has  the  best  chance  of 
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being  received  correctly.  Multi-path  fading  effects  also  cause 
random  power  variation  among  the  received  packets.  Hence,  the 
near/far  effect  and  multi-path  fading  combine  to  allow  a 
receiver  to  correctly  receive  a  packet  even  when  multiple 
packets  are  transmitted  simultaneously.  Recent  research  has 
shown  that  each  of  the  above  effects  do  improve  the  throughput 
of  an  Aloha  channel  considerably  [1-3].  However,  this  is 
achieved  in  favor  of  the  users  near  the  base  station  and  at 
the  expense  of  the  far  users.  The  far  users  have  to  retransmit 
their  packets  more  often  than  the  near  users  and  hence  suffer 
a  longer  packet  delay. 

To  correct  this  problem,  a  ring  model  in  which  the 
coverage  area  is  divided  into  3  concentric  rings  with  the  base 
station  in  the  center  is  investigated  in  this  thesis.  Users  in 
different  rings  use  different  strategies  to  retransmit  their 
collided  packets.  A  ring  that  is  far  from  the  base  station 
will  have  a  shorter  average  retransmission  delay  than  a  ring 
that  is  closer  to  the  base  station  to  compensate  for  the 
higher  capture  probability  of  the  near  users.  As  a  result, 
users  in  a  ring  model  may  achieve  approximately  the  same 
packet  delay. 

This  thesis  presents  an  analysis  of  the  near/far  effect  on 
channels  employing  variable  packet  Aloha  in  a  local/mobile 
radio  ring  model.  Unlike  idealized  Aloha  in  which  a  tagged 
packet  is  only  received  successfully  when  it  does  not  collide 
with  other  packets,  the  capture  effect  enables  the  base 
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station  to  accurately  decode  the  tagged  packet  even  when  it  is 
involved  in  a  collision.  The  result  of  the  analysis  is  used  to 
find  the  equal  riny-throughput  of  a  network  consisting  of 
three  concentric  rings. 

Chapter  II  presents  a  typical  system  model  of  local/mobil 
radio  communications.  Chapter  III  analyzes  the  dynamics  of 
number  of  interfering  packets  in  variable  packet  Aloha 
channels.  The  probability  of  capture  of  a  tagged  packet  by  the 
base  station  and  the  system  throughput  will  be  established  in 
chapter  III.  Finally,  in  chapter  IV  a  network  of  three  rings 
is  analyzed  to  find  the  ring  widths  such  that  the  throughput 
in  each  ring  is  approximate  equal. 
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II 


SYSTEM  MODEL 


The  system  model  is  a  network  consisting  of  users  assumed 
to  be  uniformly  distributed  in  an  annular  ring  surrounding  the 
base  station.  Figure  1  depicts  the  system  model.  Assume  t-hat 
the  user  traffic  density  is  g(z),  where  r  is  the  distance  to 
the  base  station.  The  density  g(r)  can  be  thought  of  as  the 
normalized  packet  traffic  per  unit  area  ac  a  distance  r  from 
the  base  station.  The  pajket  generation  rate  in  a  given  area 
depends  on  r  and  is  independent  of  direction.  Then  the  traffic 
in  a  differential  ring  of  width  dr  at  a  radius  r  from  the  base 
station  is  [4] 

G(dr)  =2nrg(x)  dr  (1) 

The  traffic  generated  within  a  distance  r  of  the  base  station 
is 

r 

G(r)  =27 ifxg(x)dx  (2) 

0 

The  total  offered  traffic  in  the  network  can  be  obtained 
as  G  =  G(r  -*  oo) 

m 

G  =  2k Jxg(x)  dx  (3) 

o 
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The  spatial  distribution  of  packet  generation  is  found  to 
be  the  probability  a  packet  is  generated  within  a  distance 
R  <  r 

Z 

FR{r)  =  PziP.zr)  =  ~  J  xg{x)  dx  (4) 

o 

Ther  fore,  the  probability  density  function  for  packet 
generation  is  found  to  be 

fR(r)  =  15 -rg(r )  (5) 

(3 


Figure  1.  The  system  model 


Assume  that  users  are  uniformly  distributed  in  an  annular 
ring  R,  <  r  i  R2  around  the  base  station,  then 

G 


g(r)  =  •(  it  (/?/-*!> 

l  0 


,  s;  r  s;  *2 
,  otherwise 


(6) 


Substituting  (6)  into  (5),  we  obtain 


fR(r) 


,  £  r  s  R2 

,  otherwise 


(7) 
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Figure  2  shows  the  distribution  of  packet  generation  of  a 
ring  with  R1  <  r  <  R2. 

g(r) 

G  _ 

tc  (Rl  -  rI) 

R1  R2  r 
Figure  2.  User  traffic  density  g(r) 

There  are  two  main  effects  that  influence  the  performance 
of  local/mobile  radio  communications,  namely  the  near/far 
effect  and  multi-path  fading.  To  model  the  near/far  effect, 
assume  that  signal  power  depends  on  the  spatial  distribution 
of  users  in  the  network.  The  mean  power  of  a  packet  at  a 
distance  r  from  the  base  station  is  of  the  general  form 

w  =  cr'“  (8) 

where  a  is  between  3  and  5  for  ultra  high  frequency  (UHF) 
applications.  The  constant  c  depends  on  transmitted  power  and 
the  gains  and  the  heights  above  the  ground  of  the  base  station 
and  mobile  antennas.  Assuming  all  users  are  identical  we  can 
normalize  the  constant  c  to  unity  for  all  users.  Hence, 

w  -  r'“  (9) 

Let  the  power  random  variable  W  be  defined  by 
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W  =  R-* 


(10) 


The  event  {W  <  w}  occurs  when  event  A  =  {R'“  <  w}  occurs. 
A  can  be  rewritten  as  A  =  {R  >  w'1/a} .  Thus, 

Fw(w)  =  P[Rzw~1/a ]  (11) 

Therefore,  (11)  is  equivalent  to 

F„(w)  =  1  -  Fr  ( W1/o )  (12) 

Differentiate  both  side  of  (12) 

f„(w)  =  -fg(w~V)  -~~a  ]  (13) 


Expanding  (13),  we  get 


-l 


-l 


fw(w)  =  fB(w  “  ) 


w 


(14) 


Substituting  (7)  into  (14) ,  we  obtain 


f  w  ( w)  = 


a  ( i?2  -  i?i ) 


w 


■t-i+i) 


(15) 


For  a  typical  value  a 
becomes 


=  4  in  UHF  communications,  (15) 


fw(w)  = 


w 


-  (3/2) 


2  (i?|-i?i) 


(16) 


One  of  the  important  factors  in  determining  the  quality  of 
signal  reception  is  the  number  of  interferers  to  the  tagged 
packet.  Thus,  it  is  of  interest  to  analyze  the  dynamics  of  the 
system.  This  is  presented  in  the  next  chapter. 
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II.  SYSTEM  DYNAMICS 


In  this  study,  we  propose  variable  packet  Aloha  instead  of 
the  conventional  fixed  length  packet  Aloha.  Variable  packets 
do  not  require  padded  bits  to  fill  up  a  packet  to  its 
specified  length  which  may  reduce  the  throughput  of  fixed 
packet  Aloha.  Variable  packets  also  do  not  require  bit 
stuffing  and  de-stuffing  logic  which  may  increase  the  cost  of 
hardware.  The  channel  traffic  rate  is  assumed  to  obey  a 
Poisson  process  with  parameter  g  packets  per  second.  The  newly 
created  packets  are  of  variable  length  x  (seconds)  and  are 
independently  and  identically  distributed  according  to  the 
probability  density  function  a(x)  with  mean  x.  All 
retransmitted  packets  have  their  length  x  redrawn  afresh  from 
a(x)  for  every  retransmission. 

The  probability  that  n  packets  interfere  (overlap)  with  an 
arbitrary  tagged  packet  during  its  transmission  is  given  by 
[5] 

Pr{n }  =  i9{x*u)J^e.g{x.u)  (17) 
n\ 

where  x  is  the  tagged  packet  length,  u  is  the  random  length  of 
a  preceding  packet,  and  n  is  the  number  of  interfering 
packets . 

The  probability  that  the  tagged  packet  encounters  no 
overlap  from  previous  packets  and  the  next  inter-arrival  is 


8 


larger  than  x  is  given  by  n  =  0.  This  is  also  the  probability 
of  successful  tagged  packet  transmission  for  the  case  where  no 
interfering  packets  can  be  tolerated  by  a  tagged  packet. 

This  analysis  considers  the  case  when  the  base  station  is 
captured  by  a  tagged  packet  in  the  presence  of  interfering 
packets  if  the  ratio  of  the  tagged  packet  signal  power  and  the 
joint  interfering  packets  signal  power  exceeds  a  threshold  yQ. 

The  number  of  interfering  packets  in  a  variable  packet 
Aloha  channel  is  a  stochastic  process  which  is  driven  by  four 
types  of  interferers.  The  first  type  of  interferers  are 
designated  as  Early-Late  packets.  These  are  packets  that  begin 
before  the  tagged  packet  and  end  after  the  tagged  packet.  The 
remaining  types  of  interferers  have  similar  descriptive 
nomenclatures  and  are  designated  as  Early-Early  packets,  Late- 
Early  packets,  and  finally  Late-Late  packets.  Figure  3  depicts 
four  types  of  interferers  relative  to  the  tagged  packet. 


TAGGED  PACKET 


EARLY -LATE  PACKET 


EARLY -EARLY  PACKET 


LATE -EARLY  PACKET 


LATE -LATE  PACKET 


Figure  3.  Four  types  of  interfering  packets 
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The  arrival  and  departure  of  n  interfering  packets 
partition  the  tagged  packet  into  several  non-overlapping 
intervals  of  random  length.  The  number  of  these  intervals, 
given  a  known  number  of  interfering  packets  n  and  their 
corresponding  type  of  interfering  packet?  found  in  is  [6] 

Total  Intervals  =  number  of  Early-Early  packets 

+  2  x  ( number  of  Late-Early  packets)  (18) 
+  number  of  Late-Late  packets  +  l 

The  relationship  in  (18)  is  based  on  the  following  facts. 
There  is  initially  one  interval,  the  tagged  packet  itself, 
when  there  are  no  interfering  packets.  The  presence  of  a 
single  Early-Early  packet  divides  the  tagged  packet  into  two 
smaller  intervals;  it  adds  one  more  to  the  total  number  of 
intervals.  The  presence  of  a  single  Late-Late  packet  also 
partitions  the  tagged  packet  and  provides  one  more  to  the 
total  number  of  intervals.  The  presence  of  a  single  Late-Early 
packet  will  further  divide  the  tagged  packet  and  adds  two  more 
to  the  total  number  of  intervals.  The  Early-Late  packet 
overlaps  entire  duration  of  the  tagged  packet  and  therefore 
contributes  no  extra  intervals  to  the  total  number  of 
intervals . 

Each  interval  has  associated  with  it  a  total  number  of 
interferers  I(ts),  where  t;  is  the  start  time  of  the  i'th 
interval  of  the  tagged  packet.  Each  unique  sequence  of  total 
interferers  per  partition  I(tj)  for  all  i  is  called  a 
realization.  Figure  4  depicts  a  specific  realization. 
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TAGGED  PACKET 


Figure  4.  An  example  of  a  specific  realization. 

The  total  number  of  inter ferers  I(tj)  =  {2 ,3, 4, 3, 2} 

For  this  specific  realization,  the  number  of  interfering 
packets  is  four.  There  is  one  of  each  type  of  packet:  one 
Early-Late  packet,  one  Early-Early  packet,  one  Late-Early 
packet,  and  one  Late-Late  packet.  The  total  number  of 
intervals  during  the  tagged  packet  according  to  (18)  is  five. 
The  total  number  of  interferers  for  each  interval  is  shown  at 
the  top  of  Figure  4.  We  are  interested  in  finding  every 
possible  realization  for  a  given  n  users.  From  those 
realizations,  the  number  of  realizations  C-  in  which  the 
maximum  number  of  interferers  present  at  any  time  during  the 
tagged  packet  interval  does  not  exceed  some  number  j,  will  be 
calculated. 

Finding  all  possible  realizations  is  very  tedious  but  is 
easily  obtained  via  computer.  TABLE  I  lists  the  algorithm  of 
the  main  program. 
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TABLE  I.  MAIN  PROGRAM  ALGORITHM 


%  N+l  is  a  number  of  users  in  the  network 
for  n  =  1:N, 

%  call  subroutine  make_comb  to  find  all  possible 
%  combination  of  interfering  packets, 
length  =  make_comb ( ) ; 
for  i  =  1: length, 

%  call  subroutine  realization  to  build  a  binary 
%  tree  of  all  possible  realizations 
realization ( ) ; 

%  call  get_cj  to  count  number  of  realizations 
%  which  has  maximum  number  of  interferer  not 
%  exceeding  j 

get_cj ( ) ; 

end 

end 


For  a  given  n  from  1  to  N,  where  N  is  the  maximum  number 
of  inter ferers  in  the  system,  subroutine  make_comb  finds  all 
possible  combinations  of  four  types  of  packets  taken  n  at  a 
time.  Basically,  subroutine  make_comb  assigns  different  values 
each  time  for  each  type  of  interfering  packets  such  that  the 
total  of  all  interfering  packets  is  equal  to  n  and  enumerates 
the  number  of  distinctive  combinations  found.  Subroutine 
make_comb  uses  the  algorithm  found  in  [7]  to  generate  all  the 
combinations.  For  example,  given  n  =  2,  the  set  of 
combinations  of  four  types  of  interfering  packets  will  be  U  = 
{(2, 0,0,0),  (0,2, 0,0),  (0,0, 2,0),  (0,0, 0,2),  (1,1, 0,0), 
(1,0, 1,0),  (1,0, 0,1),  (0,1, 1,0),  (0,1, 0,1),  (0,0, 1,1)). 

For  each  combination  of  four  packet  types  taken  n  at  a 
time,  the  subroutine  realization  is  called  to  find  the  total 
number  of  realizations.  The  root  of  the  binary  tree 
corresponds  to  the  first  interval  and  contains  the  initial 
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number  of  interferers.  In  the  first  interval,  the  interferers 
can  only  be  those  packet  types  beginning  before  the  tagged 
packet;  i.e.,  Early-Late  packet (s)  and  Early-Early  packet (s). 
The  initial  number  of  interferers  will  be  the  sum  of  those 
packets.  The  children  of  the  root  of  the  binary  tree 
corresponds  to  interval  2  and  their  children  corresponds  to 
interval  3  and  so  on  for  the  rest  of  all  intervals.  The  next 
intervals  of  the  tagged  packet  contain  either  one  packet  more 
or  one  packet  less  than  the  number  of  interfering  packets  of 
the  previous  interval  because  their  existence  depends  on 
either  the  arrival  or  the  departure  of  one  of  the  interfering 
packets.  However,  the  last  interval  should  be  equal  to  number 
of  those  packets  ending  after  the  tagged  packet  ends.  We 
choose  to  implement  the  binary  tree  by  having  the  left  child 
node  taken  care  of  the  case  of  arrival  of  an  interfering 
packet  and  the  right  child  node  for  the  case  of  departure  of 
an  interfering  packet.  In  other  words,  the  number  of 
interferers  stored  at  a  left  child  node  of  a  parent  node  is 
equal  to  the  parent's  number  of  interferers  plus  one  if  the 
parent  node  is  not  greater  than  n,  the  maximum  number  of 
interferers.  If  it  is  greater  than  n  ,  that  left  child  node  is 
set  to  NULL;  i.e.,  any  realization  containing  that  node  is 
invalid.  The  number  of  interferers  stored  at  a  right  child  is 
equal  to  the  parent  node's  minus  one  if  the  number  of 
interferers  of  the  parent  node  is  greater  than  the  minimum 
number  of  interferers  which  is  equal  to  number  of  Early-Late 
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packet  (s)  .  If  it  is  not  greater  than  the  minimum  number  of 
interferers  then  the  right  child  node  is  set  to  NULL.  Any 
realization  which  has  the  number  of  interferers  in  its  last 
interval  not  equal  to  the  sum  of  the  number  of  Early-Late 
packets  and  Late-Late  packets  is  also  invalid.  Figure  5  shows 
a  binary  tree  of  a  specific  combination  of  four  types  of 
interfering  packets  (1,1, 1,1). 


Figure  5.  A  binary  tree  implementation  of  (1,1, 1,1) 

The  subroutine  get_cj  traverses  the  binary  tree  to  collect 
all  the  valid  realizations.  The  count  of  realizations  which 
have  the  same  number  of  maximum  number  of  interferers  j  at  any 
time  is  Cj . 

TABLE  II  shows  the  output  of  the  computer  program  for  n  =1 
to  7.  For  the  case  of  n  =  1,  there  are  only  four  possible 
realizations,  namely  the  interferer  can  be  Early-Early,  Early- 
Late,  Late-Early,  or  Late-Late  packet.  Each  realization  has 
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only  one  known  interferers.  Therefore,  the  value  of  0,(1)  is 
four.  Other  values  of  Cj(l)  such  as  C2(l)  ,  C3(l)  ,  etc.  are  0. 
For  the  case  of  n  =  2 ,  there  are  10  combinations  of  four 
packet  types  taken  two  at  a  time.  These  combinations  generate 
14  unique  realizations.  Of  these  14  realizations,  four  have  a 
maximum  number  of  interferers  of  1  at  any  time,  and  10  have  a 
maximum  number  of  interferers  of  two  at  any  time.  The  value  of 
C,(2)  is  4  and  C2(2)  is  10.  Other  Cj  (2)  values  are  0.  The 
function  IT  is  the  total  number  of  all  possible  realizations 
for  n  known  interfering  packets.  The  function  f(n)  is  the 
inverse  of  IT. 


TABLE  II.  MAXIMUM  NUMBER  OF  INTERFERING  PACKETS  GIVEN  n 


n 

c? 

C4 

C<i 

CA 

c7 

.  It 

f  (n) 

n 

■1 

0 

0 

0 

0 

0 

0 

4 

0.25 

Wm 

H 

10 

0 

0 

0 

0 

0 

14 

0.0714 

1:1 

24 

20 

0 

0 

0 

0 

48 

0.0208 

Mm 

50 

75 

35 

0 

0 

0 

164 

0.0061 

KB 

■1 

100 

225 

176 

56 

0 

0 

561 

0.0018 

H 

200 

625 

664 

350 

84 

0 

1927 

0.0005 

H 

u 

400 

1698 

2250 

1565 

624 

120 

6661 

0.0002 
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III.  PROBABILITY  OF  CAPTURE  AND  THROUGHPUT  OF  VARIABLE  PACKET 


ALOHA  CHANNELS 

A.  PROBABILITY  OF  CAPTURE 

As  we  have  mentioned  earlier,  the  near/far  effect  arises 
when  users  are  scattered  at  different  distances  from  the  base 
station  and  causes  the  power  levels  of  received  packets  to 
vary  substantially.  Thus,  the  arriving  packet  with  the  highest 
power  level  has  a  best  chance  of  being  received  correctly. 
This  model  in  which  a  variable  packet  Aloha  protocol  is 
employed  has  better  throughput  than  the  ideal  Aloha  protocol 
due  to  the  near/far  effect.  A  tagged  packet  has  the  capability 
to  capture  the  base  station  even  when  it  is  involved  in  a 
collision.  We  assume  that  a  tagged  packet  will  capture  the 
base  station  if  its  signal-to  interference  ratio  exceeds  a 
threshold  yQ.  We  also  assume  that  thermal  noise  is  negligible 
since  it  is  small  compared  to  the  competing  interference  of 
interfering  packets.  In  this  section  we  calculate  the 
probability  of  capture  for  the  variable  packet  Aloha  ring 
model  and  develop  an  expression  for  the  throughput. 

Let  the  random  variable  X  denote  the  power  of  a  tagged 
packet  with  the  probability  density  function  (pdf)  fx(x),  and 
the  random  variable  Y  denote  the  joint  interfering  power  of  n 
>  1  interfering  packets  with  the  conditional  probability 
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density  function  fY(y|n).  With  the  assumption  of  independence, 
the  conditional  pdf  fy(y|n)  is  found  by  the  relationship  [6] 


fY(y\n)  =  f  (n)  ^  Cj(n)  [fw  {w)  ]  (19) 

j- 1 

where  j®  denotes  j-fold  convolution,  f(n)  is  the  inverse  of 
total  number  of  realizations  given  n  interfering  packets,  and 
Cj(n)  is  the  number  of  realizations  in  which  a  maximum  number 
of  interfering  packets  equals  j  given  n  interferers.  Both  f(n) 
and  Cj(n)  are  listed  in  TABLE  II. 

Let  the  random  variable  Z  be  the  signal-to-interference 
ratio 


Z  = 


X 

Y 


(20) 


then  the  probability  density  function  of  Z  is  [8] 

m 

fziz\n)  =  fyfX'Y(x,y\n)  dy  (21) 

0 

Given  the  fact  that  X  and  Y  are  independent  random 
variables,  we  obtain 

fz(z\n)  =  J'yfx(yz)fY(y\n)dy  (22) 

0 

The  probability  of  capture,  that  is  the  probability  that 
Z  is  greater  than  a  threshold  y0  given  n  interfering  packets 

Pr {capture\n}  =  Pr{ y0  <.  z  <>  »|n}  (23) 

is  obtained  by  integrating  f2(z|n)  from  y0  to  00 •  Hence 
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(24) 


The  overall  probability  of  capture  of  a  tagged  packet  is 
found  by  multiplying  the  conditional  probability  of  capture 
given  n  interfering  packets  by  the  probability  that  n  packets 
interfere  with  a  tagged  packet  during  its  transmission,  and 
summing  over  all  n. 

N 

Pr{captu±e }  =  ^2  Pr{capture\n}  Pr{n}  (25) 

n-0 


where  N+l  is  the  number  of  users  in  the  system. 

Substituting  (17)  into  (25),  we  obtain 

Pr  {capture}  =  V"  .  tg(^u)]  e'ar(x*u)  Pr  {capture\n}  (26) 
^  n ! 


n*0 


For  the  purpose  of  illustration,  consider  a  network  with 
R,  =  0  and  R2  =  1 .  The  pdf  of  packet  power  fu  in  (16)  becomes 

fw(w)  =  —  V(3/2)  (27) 

W  2 


Equation  (21)  and  the  -v-nditional  probability  of  capture 
F  -{capture  |  n)  with  the  assumption  of  y0  being  equal  to  6  dB  is 
obtained  numerically  using  Pro-Matlab  (see  Appendix  A  for  the 
program) .  TABLE  III  shows  the  conditional  probability  of 
capture  for  n  from  0  to  7 . 

In  the  next  section  we  will  incorporate  results  in  TABLE 
III  to  find  the  throughput  of  a  variable  packet  Aloha  channel 
in  a  ring  model. 
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TABLE  III.  THE  CONDITIONAL  PROBABILITY  OF  CAPTURE 


n 

Pr{ capture | n) 

o 

1.0000 

1 

0.2002 

2 

0.1503 

3 

0.1202 

4 

0.1018 

5 

0.0894 

6 

0.0806 

7 

0.0740 

B.  THROUGHPUT  OF  A  VARIABLE  PACKET  ALOHA  RING  MODEL 

The  channel  throughput  S,  in  packets  per  mean  packet  time 
is  [5]. 

to  oo 

*-//  (x,  u)  a  (x)  a  (u)  dxdu  (28) 

0  0 

where  a(x)  and  a(u)  are  exponential  density  functions  with 
mean  x  of  the  tagged  packet  length  and  the  preceding  packet 
length,  respectively.  The  quantity  0(x,u)/2  is  the  conditional 
channel  departure  rate  [5], 

-i-4>  (x,  u)  =  -igr(x  +  u)  Pr  (capture)  (29) 

Substituting  (26)  and  (29)  into  (28),  we  have 
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s-//-|g(**u)  ,30, 

0  0  n*° 

.  Pr  {capture |n}  a  (x)  a  ( u )  dxdu 

Utilizing  the  binomial  identity 

( x  +  u)n  =  ^  I  ?)xk  un~k  (31) 

kTo  \  kl 

and  expanding  (30) ,  we  obtain 

M  « 

S  ~  —f  f  [gxe'9*  e_?u  +  gue'911  e_?*  ]  ^  Pr{capture|.n} 

2  Vo  D-°  ■n!  (32) 

( u}x k  un~k  a a ( ^^u 

k-o  \ 

Further  manipulating  terms  in  the  above  equation  gives 

•o 

S  =  -i-g]T  "^T  Pr (capture I72) 

^  /)-0  ^  * 

J3  *  “ 

•  |  frjf  xk*1e'gxa(x)  dx.  j  un~ke~gua(u)  du  (33) 

U  *  "* 

+  £(?)[  u^e-^a(u)du.  f  xke~9Xa{x)  dx 

k-o\klJ0  o 

We  assume  that  the  packet  length  is  exponentially 
distributed  and  the  average  length  is  x.  Then  we  have 


a  (x) 

=  -Le~x/x 

(34) 

X 

a  (u) 

=  - Le’u/X 

X 

(35) 

Define  the  average  traffic  rate  G  as 

G  =  gx 

We  apply  (34),  (35),  (36)  and  (9]. 
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j  xn  e~  **x  dx  =  73!ji'(n+1> 


to  solve  (33) 


The  first  integral  term  is  evaluated  as 


J  x<,k  +  l)  e-?xa  (x)  dx  = 


If  v(Jc*l>0-(y+T)x, 


U£±±LL  (g+i)-U+2> 

.  L£ll2  i  (G+l)-<*+2> 


Similarly,  the  second  integral  is  found  to  be 


J  u{n'k)  e'^aiu)  du  = 


=  -kf  u{n~k)  e  tsr*"ff)udu 
uJ 

=  (n-_k)  !  (  -(n-*+1) 

x  y  x 

=  .{JzlFL!-  ((?+i) -<»-**!> 

y  ~~  i.  n  Jc ) 


The  third  integral  term  is 


?u(n-*+1)  e_sru  a(u)  du  =  k  +k11).i-  (G+ 1) 

J  ZT~{n-k*X) 


The  final  integral  term  is  reduced  to 


f  xne-**a{x)dx  =  —  (  G+l) -(jc+1> 

J  x-k 


Combining  all  integral  terms  in  (38),  (39),  (40),  and  (41) 


into  (33) ,  we  get 


s  -  j^Prlcapture\n}  ■(G^,).'>1°0’  £  n!  (nt2 

^  n-o  A  Jc-o 


)  (42) 
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Applying  (36)  and  the  equality 

n\ (n+2)  =  n\  (n+l)  (n+2)  (43) 

Jc  ■  0 

and  further  simplifying  (42) ,  we  obtain 

S  =  Y  +  — -  Prlcapture\n}  G"*1  (G+l)  ~<n+3)  (44) 

n«0  2 

For  the  network  with  R1  =  0,  Rj  =  1,  and  the  conditional 
probability  of  capture  in  Table  III,  throughput  is  shown  in 
Figure  6. 


0.25 


Figure  6.  Throughput  of  a  ring  model 
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IV.  APPROXIMATE  EQUAL  RING  THROUGHPUT 


In  ♦■.his  section,  re  expand  cur  analysis  to  a  system  of 
multiple  rings.  As  we  mentioned  before,  the  probability  of 
capture  favors  near  users  at  the  expense  of  far  users 
resulting  in  lower  throughput  for  the  far  users.  A  multiple 
ring  system  compensates  for  this  deficiency  by  using  separate 
retransmission  strategies  for  users  in  different  rings.  A  user 
in  a  ring  that  is  far  from  the  base  station  is  given  a  shorter 
average  transmission  delay  than  a  user  in  a  ring  closer  to  the 
base  station  to  compensate  for  a  higher  capture  probability  of 
the  near  ring. 

A.  PROCEDURE  TO  FIND  EQUAL  RING  THROUGHPUT 

A  procedure  analogous  to  the  procedure  found  in  [3]  is 
used  to  find  the  ring  width  for  each  ring  such  that  the 
throughput  in  each  ring  is  approximately  equal.  The  procedure 
consists  of  the  following  iterative  process: 

1.  Select  initial  ring  widths  Rf  =  (r,.  <  r  <  ri+1,},  where 
i  =  1,2, ...,M. 

2.  Find  the  conditional  ring  capture  probability 
Pr{ capture | n,Rs } (see  Appendix  B) . 

3.  Find  the  traffic  rate  Gf  in  each  ring 

Gi  =  G[FR(rUl) -FR(rd)]  (45) 
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where  i  =  and  FR(r,.)  is  given  in  (B3) 

4 .  Find  the  ring  throughput  for  each  ring  from  (Appendix 

C)  . 

=  ^-ri{captus:ii\n,Ri)GiGQ(G+l)-u*y,\  (46) 

nTo  L  2  J 

5.  If  S,  «  S2  «. . .«  SM  for  0  <  S{  <  min{max  Sj}  go  to  next 
step.  Otherwise,  go  back  to  step  1.  For  our  computational 
purpose,  the  maximum  relative  difference  of  ±0.05  is 
acceptable  for  the  approximate  equality. 

6.  Calculate  the  total  throughput  of  the  model  S  =  SSj. 
For  illustrative  purpose,  consider  a  system  of  three  rings 

R,  =  {0  <  r  <  0.52},  =  {0.52  <  r  <  0.8}  and  R3  =  {0.8  <  r  < 

1}.  The  threshold  y0  is  taken  to  be  6  dB.  Fig  7  shows  the  ring 
throughput  of  ring  R1#  R^  Rj  and  the  total  ring  throughput  of 
the  model.  The  ring  throughputs  are  approximately  equal  for  0 
<  G  <  0.25. 

B.  AVERAGE  RETRANSMISSION  DELAY 

To  determine  the  average  retransmission  delay  Tj  for  ring 
Rj,  we  need  first  to  define  a  retransmission  strategy  when  a 
collision  occurs.  One  suggestion  [10]  is  to  choose  an 
arbitrary  time  interval  and  select  a  uniformly  distributed 
random  retransmission  time  within  that  interval.  Let  the  time 
interval  covers  Kj  message-unit  times  where  each  message-unit 
is  an  average  x  units  of  time.  Then  retransmission  takes  place 
within  1  to  K(  time  units  after  a  collision  has  occurred.  The 
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Figure  7.  Approximate  Equal  Ring  Throughput 
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approximate  number  of  retransmissions  is  found  to  be 


1  -PxicaptuielRj) 
PrTcapture\R 


where  Pr { capture  j  R( }  is  the  probability  of  capture  in  ring  R( 
per  mean  packet  time  is  found  by  (Appendix  D) 


Pi{capture\Rj)  =  ^  Pr{capture\n, Rj  (n  +  1)  Gn  (G+l)  ~in*2)  (48) 

12-0 

The  average  retransmission  time,  ignoring  propagation 
delay,  is  approximately 


Tt  «  O.BEUCi  +  l) 


By  fixing  the  value  K„  of  the  outermost  ring,  the 
remaining  value  K.  can  be  obtained  by  letting  T;  =  TM  and 
substituting  (47)  into  T{  and  TM  expressions,  we  obtain 


K  _  Px{capture\Ri}  {l-Pi{captuie\RM})  (km+i) 

1  Pr{capture\RM)  (l  -Pr {capture ji^} )  x  x  j 


TABLE  IV  shows  the  number  of  retransmission  K1  and  K2  of 
ring  R1  and  ring  R2,  given  Kj  =  8  of  x  time-units.  For  o  <  G  < 
Gmax  =  °-25/  ^e  average  K,  is  29.37  and  Kj  is  8.70. 
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TABLE  IV.  VALUES  OF  K, 


G 

Ki 

K? 

0.05 

27.9389 

8.6379 

8 

0.10 

28.6474 

8.6694 

8 

0.15 

29.3634 

8.7003 

8 

0.20 

30.0863 

8.7306 

8 

0.25 

30.8158 

8.7603 

8 
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V.  CONCLUSION 


We  have  evaluated  the  probability  of  capture  taking  into 
account  the  near/ far  effect.  These  results  are  then  used  to 
compute  the  throughput  of  a  ring  model  utilizing  the  variable 
Aloha  protocol  in  a  local/mobil  radio  communications.  We  have 
investigated  a  ring  model  in  which  different  retransmission 
strategies  are  used  to  alleviate  the  near/far  effect  problem. 
By  assigning  shorter  average  retransmission  delay  to  users 
farther  from  the  base  sitation  and  longer  average  transmission 
delay  to  users  nearer  to  the  base  station,  all  users  may 
achieve  approximately  the  same  average  retransmission  delay. 
In  this  thesis,  we  only  investigated  a  system  of  three  rings. 
Further  work  is  needed  to  determine  the  number  of  rings  that 
produces  the  maximum  throughput.  The  procedure  for  finding 
equal  ring  throughput  is  a  process  of  trial  and  error  and 
requires  greater  effort  when  a  number  of  rings  is  greater  than 
three.  The  capture  measure  using  signal-to-interference 
threshold  provides  a  less  accurate  prediction  of  channel 
throughput  than  does  a  capture  measure  using  packet  correction 
reception  probability  [1],  However,  employing  packet 
correction  reception  probability  as  a  capture  measure  for 
variable  packets  is  very  difficult.  It  is  an  unwielding  task 
to  find  the  probability  of  bit  error  for  variable  length 
packets. 
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APPENDIX  A.  COMPUTER  PROGRAM  LISTINGS 


A.  PROGRAM  TO  FIND  THE  MAXIMUM  NUMBER  OF  INTERFERING  PACKETS 

/*  The  program  is  written  in  Unix  C  to  find  the  number  */ 
/*  of  realizations  Cj  which  has  the  maximum  number  of  */ 
/*  interferers  equal  to  j .  It  first  calls  subroutine  */ 
/*  make_comb  to  find  all  the  combinations  of  four  types  */ 
/*  of  interfering  packets  taken  n  at  a  time.  Then,  for  */ 
/*  each  combination  it  calls  subroutine  realization  to  */ 
/*  build  a  binary  tree  of  realizations.  */ 

/*  Finally,  subroutine  get_cj  is  called  to  find  Cj .  */ 

#include  <stdio.h> 

#define  comb_lgth  230 
#define  type  4 
#define  cj_lgth  10 
#def ine  INVALID  -1 

#define  NODE_SIZE  sizeof (struct  node) 

/*  Global  variables  declaration  */ 

struct  node 
{ 

int  num_intfr;  /*  number  of  interfering 

packets  */ 

struct  node  *left;  /*  left  pointer  of  a  node  in  a 

binary  tree  */ 

struct  node  *right;  /*  right  pointer  of  a  node  in 

a  binary  tree*/ 

/*  file  descriptor  */ 

/*  number  of  users  in  system  */ 

/*  number  of  intefering  packets  in 
the  first  interval  i.e.  number  of 
early  beginning  packets  */ 

/*  number  of  late  ending  packets  */ 
/*  number  of  maximum  possible 
interfering  packets  i.e.  number  of 
users  */ 

/*  minimum  number  of  interfering 
packets  */ 

/*i.e.  number  of  early-late  packets  */ 
int  num_interval ;  /*  number  of  intervals  of  tagged 

packet  */ 

int  cj [cj_lgth] [cj_lgth] ;  /*  array  of  Cj  */ 


FILE  *file; 
int  num_user; 
int  init  intfr; 


int  final_intfr; 
int  max  intfr; 


int  min  intfr; 
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/*  Function  declaration  */ 
int  make_comb ( ) ; 
struct  node  *realization() ; 
void  get_cj(); 
void  print_out() ; 
char  *malloc() ; 


main() 


{ 


int  i ,  j  ; 

int  length;  /*  length  of  array  comb  */ 

int  comb [comb_lgth] [type] ;/*array  of  combination  of 

possible  packet  types*/ 


struct  node  *root; 


file  =  fopenC'datafile”, "a")  ; 

root  =  (struct  node  *)malloc (NODE_SIZE) ; 

root~>left  =  NULL; 

root->right  =  NULL; 

/*  clear  array  cj  */ 
for ( i=0 ; i<c j_lgth ; i++) 

{ 

for ( j=0; j<cj_lgth; j++) 
cj [i] [j]  =  0; 

) 

for  (num_user=l  ;num__user<=10  ;num_user++) 

{ 

length  =  make_comb (comb) ; 
for(i=0;  i<=length;  i++) 

{ 

init_intf r  =  comb [ i ] [ 0 ] +comb [ i ] [ 1 ] ; 
final_intfr  =  comb[ij [0]+comb[i] [type-1] ; 
max_intfr  =  num_user; 
min_intfr  =  comb[i][0]; 
num_interval  =  comb[i][l]  + 

2*comb[i] [2]+comb[i] [type-1 ]+l; 
root->num_intfr  =  init_intfr; 
root  =  realization (root) ; 
get_cj (root) ; 

) 

) 

print_out() ; 
fclose(file) ; 


/*  The  subroutine  make_comb  is  to  find  all  possible  */ 
/*  combinations  of  four  types  of  interfering  packets  */ 
/*  taken  n  at  a  time.  */ 


int  make_comb (comb) 

int  comb [comb_lgth] [type] ; 


n 


int  pkt_tvpe [ type ] ; 
int  length  =  0? 
int  index_track  =  -1; 
int  num_track; 
int  i,  j; 

pkt_type[0]  =  nmn_user; 
num_track  =  num_user ; 

comb[0][0]  =  pkt_type[0];  /*  first  combination  */ 

for  (i=l;  i<type;  i++) 

pkt_type[i]  =  comb [ length] [i]  =  0; 
while (pkt_type[type-l]  !=  num_user) 

{ 

if(num_track  >1) 
index_track  =  -1; 
index_track++  ? 

num_track  =  pkt_type[index_track] ; 
pkt_type[ index_track]  =  0; 
pkt_type[0]  =  num_track  -  1; 
pkt_type  [  index_tmack  +  1  ]  ++ ; 
length++ ; 

for ( j=0 ;  j<type;  j++) 

comb [ length ] [ j ]  =  pkt_type [ j ] ; 

) 

return (length) ; 

} 

/*  Subroutine  realization  employs  subroutine  build_tree  */ 
/*  to  build  a  binary  tree  of  realizations.  */ 

struct  node  *realization (ptr) 
struct  node  *ptr; 

{ 

int  interval; 
int  tree_level; 
struct  node  *build_tree ( ) ; 
void  build_node ( ) ; 

interval  =  1; 

tree_level  =  1; 

while ( interval<num_interval ) 

{ 

ptr  =  build_tree (ptr, interval, tree_level) ; 
interval++ ; 

) 

return (ptr) ; 
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/••  Subroutine  build_tree  searches  for  the  valid  paths  of 
realizations  */ 

struct  node  *bu*’ ld_tree (ptr ,  interval,  tree_level) 
struct  node  *ptr  ; 
int  interval? 
int  tree_level; 

{ 

if (tree_level<interval) 

{ 

if ( (ptr->left) ->num_intfr  !=  INVALID) 

{ 

tree_level++ ; 

ptr->left 
build_tree fptr->left , interval, tree_level) ; 
tree_level — ; 

) 

if ( (ptr->right) ->num_intfr  !=  INVALID) 

{ 

tree_level++ ; 

ptr->right 
build_tree (ptr->right , interval , tree_level) ; 
tree_level — ; 

} 

} 

else 

build_node (ptr, interval) ; 

return (ptr) ; 

} 

/*  Subroutine  build_node  is  to  fill  nodes  of  a  binary  */ 
/*  tree  with  the  number  of  interferers.  A  left  child  */ 
/*  node  has  the  number  of  interferers  equal  to  the  */ 

/*  parent's  plus  ne  if  the  parent's  is  not  greater  */ 

/*  than  n,  the  maximum  number  of  allowed  interferers.  */ 
/*  If  the  parent  node  is  greater  than  n,  the  left  */ 

/*  child  node  is  set  to  NULL,  i.e.  invalid.  A  right  */ 

/*  child  node  will  be  equal  to  the  parent's  number  of  */ 
/*  interferers  minus  one  if  the  parent's  is  greater  */ 

/*  than  the  minimum  number  of  interferers,  which  is  */ 

/*  equal  to  number  of  Early-Late  packets.  */ 

/*  If  it  is  not,  the  right  child  is  */ 

/*  set  to  NULL.  The  final  interval  should  have  the  */ 

/*  number  of  interferes  equal  to  the  sum  of  number  */ 

/*  of  Early_Late  and  Late-Late  packets.  */ 

/*  If  it  is  not,  that  realization  is  invalid.  */ 


void  build_n^de (ptr, interval, 
struct  node  *ptr; 
int  interval; 
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{ 


if (ptr->left  ==  NULL) 

{ 

if ( (ptr->left  =  (struct  node  *)  malloc  (NODE_SIZE) ) 
==  NULL) 

{ 

print f ("ERR0R1! ! !") ; 
exit(l) ; 

} 

(ptr->left) ->left  =  NULL; 

(ptr->left) ->right  =  NULL; 

} 

if (ptr->right  ==  NULL) 

{ 

if ( (ptr->right  =  (struct  node  *)malloc (NODE_SIZE) ) ==0) 

{ 

printf ( "ERROR2 !!!"); 
exit (1) ; 

} 

(ptr->right) ->left  =  NULL; 

(ptr->right) ->right  =  NULL; 

} 

/*  install  left  child  node  */ 
if (ptr->num_intfr<max_intfr) 

{ 

(ptr->left) ->num_intfr  =  ptr->num_intfr  +  1; 

/*  invalid  final  interferers  */ 

if  ( ( interval=num_interval-l)  &&  ( (ptr->left)  ->num_intf r 
!=f inal_intfr) ) 

(ptr->left) ->num_intfr  =  INVALID; 

} 

else 

(ptr->left) ->num_intfr  =  INVALID; 

/*  install  right  child  node  */ 
if (ptr->num_intf r>min_intfr) 

{ 

(ptr->right) ->num_intfr  =  ptr->num_intfr  -  1; 

/*  invalid  final  interferers  */ 

if ( ( interval==num_interval-l) && ( (ptr->right) ->num_intfr ! 
=f inal_intfr) ) 

(ptr->right) ->nuin_intfr  =  INVALID; 

) 

else 

(ptr->right) ->num_intfr  =  INVALID; 


/*  Subroutine  get_cj  is  to  count  the  number  of  */ 
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/*  realizations  which  has  the  maximum  number  of  */ 
/*  interferes  not  greater  than  j .  */ 

void  get_cj (ptr) 
struct  node  *ptr; 

{ 

int  tree_level; 
int  max; 

void  count_c j ( ) ; 

tree_level  =  1; 

max  =  ptr->num_intfr ; 

count_cj (tree_level,max,ptr) ; 

) 

void  count_cj (tree_level ,max, ptr) 
int  tree_level; 
int  max; 

struct  node  *ptr; 

{ 

int  old_max; 

if (tree_level<num_interval) 

{ 

if ( (ptr->left) ->num_intfr  1=  INVALID) 

{ 

tree_level++ ; 
old_max  =  max; 

if ( (ptr->left) ->num_intfr>max) 
max  =  (ptr->left) ->num_intfr ; 
count_cj (tree_level,max,ptr->left) ; 
tree_level — ; 
max  =  old_max; 

) 

if ( (ptr->right) ->num_intfr  !=  INVALID) 

{ 

tree_level++ ; 

count_cj (tree_level,max,ptr->right) ; 
tree_level — ; 

) 

} 

else 

cj [max_intfr-l] [max-1 ]++; 
return ; 

} 

void  print_out() 

{ 

int  total_rlzt; 
double  fn; 
int  i , j ; 

for ( i=0 ; i<cj_lgth ; i++) 
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{ 

fprintf (file, "%d\n" , i)  ; 
total_rlzt  =  0; 
for ( j=0; j<cj_lgth; j++) 

( 

fprintf (file, M%d  M  ,cj[i][j]); 
total_rlzt  +=  cj(i][j]; 

} 

fn  =  1.0/ (double) total_rlzt; 

fprintf (file, "%d  %8e\nM , total_rlzt, fn) ; 

} 

} 


B.  PRO-MATLAB  M-FILE  FOR  THROUGHPUT  OF  A  RING  MODEL 

The  following  M-file  is  written  to  run  on  Pro-Matlab  to 

find  the  probability  of  capture  and  throughput  of  a  ring  model 

with  R,  =  0  and  R2  =  1. 

%  set  w  from  1  to  infinitive 
w  =  1:0.1:52.2; 
resolution  =  0.1; 
range  =  1/2; 

%  find  the  pdf  of  power  of  a  packet 
fw  =  range. * (w. A (-3/2) ) ; 

%  find  j-fold  convolution  of  interfering  packets 
CONV ( 1 , : )  =  fw; 

CONVFW  =  f ft ( fw) ; 
for  i=  2:10, 

CONV ( i , : )  =  resolution  *ifft(fft (CONV(i-l, : ) )  .*  C0NVPW) ; 
end 

%  get  f(n)  and  Cj  from  data  file 
data 

for  i  =  1:10, 
for  j  =  1:513, 
tmp(j)  =  0; 
end 

for  k  =  l:i, 

tmp  =  tmp  +  (cj(i,k)  .*  C0NV(k,:)); 
end 

f y ( i , : )  =  fn(i) . *  tmp; 
end 

%  find  fz  in  fz(z/n)  =  fz*zA(-3/2) 
for  i  =  1:10, 

sum  =  range  *  (w(l) A (-0. 5) )  *  fy(i,l); 
for  j  =  2:2:512, 

sum  =  sum  +  4*range* (w( j ) A (-0. 5) ) *fy ( i , j ) ; 
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end 

for  j  =  3:2:511, 

sum  =  sum  +  2*range* (w( j ) * (-0. 5) ) *fy (i, j ) ; 
end 

sum  =  sum  +  range* (w(513) A (-0. 5) ) *fy (i, 513} ; 

fz(i)  =  (resolution/3) *sum; 

end 

%  find  the  probability  of  a  tagged  packet 
pr  =  2  *  (6A (-1/2) )  . *  fz; 

%  Throughput 
g  =  0 : 0 . 1 : 5 ; 
for  i  =  1:51, 
s ( i )  =  g(i)/((g(i)+l)A3) ; 
for  n  =  1:10, 

s(i)  =  s(i)  +  (n+1) * (n+2)/2*pr (n)  *(g(i) A (n+1) ) 
*  ( (g(i)+l) A (-(n+3) ) ) ; 

end 

end 


C.  PRO-MATLAB  M-FILE  TO  FIND  EQUAL  RING  THROUGHPUT 

This  M-file  is  written  for  Pro-Matlab  to  find  the  equal 

ring  throughputs  in  chapter  4 . 

%  set  w  from  1  to  infinitive 
w  =  1:0.05:26.6; 
resolution  =  0.05; 
range (4)  =  1/2; 

%  pdf  of  power  of  an  interfering  packet 
fw  =  1/2 . * (w. A (-3/2 ) ) ; 

%  calculate  the  j-fold  convolution  of  interfering 

%  packets 

CONV ( 1 , : )  =  fw; 

CONVFW  =  fft (fw4) ; 
for  i=  2:10, 

C0NV4 ( i , : )  =  resolution.*  if ft (fft (C0NV4 (i-1, : ) )  .* 

CONVFW ) ; 

end 

%  get  f(n)  and  Cj  from  data  file 
data 

%  find  fy 
for  i  =  1:10, 

for  j  =  1:513, 
tmp(j)  =  0; 

end 

for  k  =  l:i, 

tmp  =  tmp  +  (cj(i,k)  .*  C0NV4(k,:)); 
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end 

fy(i, :)  =  fn(i).*  tmp; 
end 

%  polynomial  fit  for  fy(n) 
for  i  =  1:10, 

fit (i, : )  =  polyfit(w, fy(i, : )  ,9)  ; 
end 

%  set  range  of  ring  1  from  0  to  0.52 
rl (1)  =  0; 
r2 (1)  =  0.52; 

range (1)  =  l/(2*(0.52*2) )  ; 

%  set  range  of  ring  2  from  0.52  to  0.8 
rl (2)  =  0.52; 
r2 (2)  =  0.8; 

range ( 2 )  =  1/(2*(0.8A2  -0.52A2)); 

%  set  range  of  ring  3  from  0.8  to  1 
rl (3)  =  0.8; 
r2 ( 3 )  =  1 ; 

range ( 3 )  =  1/(2*(1  -  0.8A2)); 

%  find  fz  in  fz(z/n)  =  fz*zA(-3/2) 
for  ring  =  1:3, 
for  i  =  1:10, 
for  z  =  6: . 1 : 24 . 4 , 
if  ring  ==  3, 
start  =  1; 
else 

start  =  r2 (ring) A (-4)/z ; 
end 

if  ring  ==1, 
ends  =  26.6; 
else 

ends  =  rl (ring) A (-4)/z ; 
end 

temp  =  0; 
if  start  <  1, 
start  =  1; 
end 

if  ends  >  26.6, 
ends  =  26.6; 
end 

if  ends  >  1, 

step  =  (ends  -  start) /8; 
y  =  start: step: ends; 
clear  f; 

f  =  polyval (fit (i, : ) ,y) ; 

sum  =  range(ring)  *  (y (1) A (-1/2) )  *  f(l); 
for  j  =  2:2:8, 

sum  =  sum  +  4*range (ring) * (y ( j ) A (-0. 5) ) *f ( j ) ; 
end 

for  j  =  3:2:7, 

sum  =  sum  +  2*range(ring) * (y ( j ) A (-0. 5) ) *f ( j )  ; 
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end 

sum  =  sum  +  range(ring) * (y (9) A (-0 . 5) ) *f  (9)  ; 

temp  =  (step/3) *sum; 

end 

if  z  ==  6, 

pr(ring,i)  =  temp  *  zA(-3/2); 
elseif  z  <  13.7, 

pr(ring,i)  =  pr(ring,i)+  2*temp*zA (-3/2 )  ; 
elseif  z  ==  13.7, 

pr(ring,i)  =  pr (ring, i)+temp*zA (-3/2)  ; 

end 

end 

pr(ring,i)  =  pr (ring, i) *. 1/2 ; 
end 

%  calculate  p(3,i)  for  z  =  13.7  to  infinitive, 
if  ring  ==  l, 
for  i  =  1:10, 
y  =  1:. 1:26. 6; 
clear  f; 

f  =  polyval ( f it ( i , : ) , y )  ; 

sum  =  range(ring)  *  (y (1) A (-1/2) )  *  f(l); 
for  j  =  2:2:256, 

sum  =  sum  +  4*range(ring) * (y ( j ) A (-0 . 5) ) *f ( j ) ; 
end 

for  j  =  3:2:255, 

sum  =  sum  +  2*range(ring) * (y ( j ) A (-0. 5) ) *f ( j ) ; 
end 

sum  =  sum  +  range (ring) * (y (257) A (-0 . 5) ) *f (257) ; 

temp  =  ( . 1/ 3 ) *sum ; 

temp  =  2* (13 . 7A (-1/2) ) *temp; 

pr(ring,i)  =  pr (ring, i) +temp; 

end 

end 

%  Throughput 
g  =  0:0.05:4; 
for  i  =  1:81, 

s(ring,i)  =  1/ (2*range (ring) ) *g (i) / ( (g (i) +1) A3) ; 
for  n  =  1:10, 

s(ring,i)  =  s(ring,i)  +  (n+1) * (n+2)/2*pr (ring, n) 

/ (2*range (ring) )*g(i)*(g(i)A(n))*((g(i)+l)A(-(n+3)) ) 
end 
end 
end 
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APPENDIX  B.  COMPUTATION  OP  CONDITIONAL  CAPTURE  PROBABILITY 
IN  RING  Rj 


Let  us  consider  a  tagged  packet  in  ring  R;  bounded  by  two 
concentric  circle  of  radii  r(  and  ri+1,  respectively.  The 
conditional  probability  density  function  for  packet  generation 
in  ring  Rf  is  given  by 


Fr  ~Fr 

0 


,  ri  i  r  sri+1 
,  otherwise 


(Bi) 


The  probability  a  packet  is  generated  within  a  distance 
r  <  r.  is 


FR(ri) 


x - ® — ~dx 

it  ( Ri-R ?) 


( B2 ) 


Solving  (B2) ,  we  have 

FrU,) 


( B3 ) 


Similarly, 


Fr  ( ri  ♦  i )  = 


rI-rI 


( B4 ) 


Substituting  (B3)  and  (B4)  into  (Bl) ,  we  obtain  the  pdf  of 
packet  generation  in  ring  R;  as 
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=  \  r^-rj 


2  r 


zi  s  rs  ri+1 


(B5) 


0  ,  otherwise 

The  corresponding  conditional  pdf  for  packet  power  X  in 
ring  Rj  can  be  obtained  from  (16)  as 


f*(*!*i)  = 


-<4> 


i  sri4l 


2  (r|+1  -r\) 

0  ,  otherwise 


(B6) 


The  pdf  of  the  random  power  of  n  interfering  packets 
anywhere  in  the  network  is  obtained  from  (19) .  Therefore,  the 
conditional  pdf  of  signal  to  interference  ratio  in  ring  Rj  is 

fz{z\n,Ri)  =  Jyfx(yz\Ri)  fY[y\ri)  dy  (B7) 

o 


Then  the  conditional  capture  probability  in  ring  R1  is 

M 

Pr{capture\n,Ri }  =  jf z  (zln,^)  dz 

To 
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APPENDIX  C.  DERIVATION  OP  THROUGHPUT  OP  A  RING  IN  A  SYSTEM 
OF  MULTIPLE  RINGS. 


The  derivation  of  the  ring  throughput  for  a  multiple  ring 
system  is  similar  to  the  derivation  in  section  II. B. 

The  channel  throughput  Sj,  in  packets  per  mean  packet  time 
becomes  [ 5 ] 

oa  am 

w/  —  $i(x,  u)  a(x)  a(u)  dxdu  (Cl) 

0  0 

where  a(x)  and  a(u)  are  exponential  density  functions  with 
mean  x  of  the  tagged  packet  length  and  the  preceding  packet 
length.  The  quantity  0j(x,u)/2  is  the  conditional  channel 
departure  rate  in  ring  R.  [5]. 

~^i(x,u)  =  -|-gri(x+u)  Pi  (capture)  (C2) 


where  gs  is  a  parameter  in  packets  per  second  of  a  Poisson 
process  of  the  channel  rate  in  ring  Rf . 

Substituting  (26)  and  (C2)  into  (Cl) ,  we  have 


Sj  £  .-*»->  [y<Y,u)r 

0  0  ^  n"0  n ' 


( C3 ) 


.  Pr{capture\n]  a  (x)  a(u)  dxdu 


Utilizing  the  binomial  identity 


(x+u) n 


( C4 ) 
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and  expanding  (C3) ,  we  obtain 

Si  =  —  pJ*  f  [xe'9*  e_ffu  +  ue e'9*  ]£}  Pr{capture\n } 

2  o  0  n-°  n‘  (C5) 

'53  (  un~k  a  (x)  a  (u)  dxdu 

kT'o  \  I 

Further  manipulating  terms  in  the  above  equation,  we  get 

m 

si  =  4^53  44  Pr  {cap ture  I73} 

^  /2-0  ^  • 

12  m  m 

•  '£[l)fxk*1e-gxa(x)dx.j  un-ke-^a(u)du  (C6) 
ic»o  *  1  q  0 

^  OO  «S 

+  53  (  u)/  un_;t+1  e'^a  (u)  du.  J  x ke~gxa  (x)  dx 

We  have  assumed  that  the  packet  length  is  exponentially 
distributed  and  the  average  length  is  x.  Hence 


a  (x)  =  -i.e'x/x 

(C7 ) 

X 

a(u)  =  -ke-u/x 

X 

( C8 ) 

Define  the  average  traffic  rate  G  of  the  system  as 

G  =  gx  ( C9 ) 

We  apply  (C7) ,  (C8) ,  (C9)  and  the  following  equation  [9]. 

j x n  e"  **x  dx  =  n !  n' <n+1>  (CIO) 


to  solve  (C6) 

The  first  integral  term  can  be  evaluated  as 
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i)  e  <sr+T )x<±x 


I  X(W)  e-^a(x)dx  =  x<*+ 

=  -ilxiU  (a+i)'(^2> 

=  (G+l) 

Similarly,  the  second  integral  is  found  to  be 


J  ula~k)  e'^aiu)  du  = 


if  u<-*>  e-|s,^"dl. 
xj 

=  <”-_*>  !  (g+i)'<n~*+1) 

X  X 

-  (n-k)  !  / *  \  - (n-k*i) 

X-(n-Jc)  V  ' 


The  third  integral  term  is 


ju(a-k*v  e-9u  a  ( u)  du  =  ■Jct1\1  (G+l) ’(n'Jc*2) 

J  v.-(n-Jc*l) 


The  final  integral  term  is  reduced  to 


/ 


xne~9Xa  (x)  dx  =  ( G+l) 

x'k 


(k+ 1) 


Combining  all  integral  terms  in  (Cll) ,  (C12) ,  (C13) 

( C14 )  into  (C6) ,  we  obtain 

si  *  -J0iE  -gf  fr(e4pture|n>  £  n!  (n*2) 

^  n-o  -X  Jc-o 

Applying  (C9) ,  the  equality 

n! (n+2)  =  n\ (n+1) (n+2) 


k-o 


(Cll) 

(C12 ) 

(C13) 

(C14) 

,  and 

(C15) 

(C16 ) 
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nd  the  fact  that  =  gfx,  and  further  simplifying  (Cl 5) ,  we 
obtain 

Si  =  pi{capture\ri)  GiGn(G  +  i) ' (n0)  (C17) 

n-o  2 


45 


APPENDIX  D.  DERIVATION  OF  PROBABILITY  OF  CAPTURE. 


The  probability  of  capture  in  a  ring  Pr{ capture},  per  mean 
packet  time  can  be  found  by 


Pr(capture)  =  f  f  £  Pr  {capturein}  Pr  [n] 

0  o  n*° 

,  a  (x)  a  ( u)  dxdu 


(Dl) 


where  a(x)  and  a(u)  are  exponential  density  functions  with 
mean  x  of  the  tagged  packet  length  and  che  preceding  packet 
length. 

Substituting  (26)  into  (Dl) ,  we  have 

90  90  mo 

Pr{capture }  =  f  f£e~fflx*u)  UU.., ,  (D2) 

0  on*o  n! 

.  Pr[capture\n]  a(x)  a  (u)  dxdu 

Utilizing  the  binomial  identity 

(x+u)”  =  un-k  (D3) 

Jt-0  \  k  I 

and  expanding  (D2) ,  we  obtain 

«  « 

Pr  [capture]  =  f  f  -2-yPr  [capture\n]  e'£  [x*u) 

o  0n-°  n]  ( D4 ) 

•53  (  ^)x  *  u”'*  a  (x)  a  (u)  dxdu 

jfro  \  kl 

Further  manipulating  terms  in  the  above  equation,  we  get 
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Pi  {capture}  =  -2—  Pi  {capture\n} 


f  xke~9*a{x)  dx.  j  un~k e'9” a  (u)  du 

JC"0  '  '  n  n 


We  have  assumed  that  the  packet  length  is  exponentially 
distributed  and  the  average  length  is  x.  Hence 


a(x)  =  -ie'x/x 

(D6) 

X 

a  (u)  =  ie*u/x 

(D7) 

Define  the  average  traffic  rate  G  of  the  system  as 

G  =  gx 

We  apply  (D6) ,  (D7) ,  (D8)  and  the  following  [9] 

jxn  e~  **  dx  - 
o 

to  solve  (D5) 

The  first  integral  term  can  be  evaluated  as 

m  m 

j  x  ke~gxa  (x)  dx  =  -L j  xke  i9*~*‘xdx 

o  X  0 

,  !  ,  1  > 

=  k-=,  (g  + 
x  x 

=  k-J~  (G  +  l)  “  tJc*1) 
xk 

Similarly,  the  second  integral  is  found  to  be 

«• 

f  uin-k)  e'^a(u)  du  =  1  (g+l) 

J  y-la-k) 


(DIO) 


(Dll) 
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Combining  all  integral  terms  in  (DIO) ,  (Dll) ,  into  (D5) , 


we  get 


“  rr  n  ( f  .  -i  \  -  in*2)  JL 

Pr{capture }  =  Pi {captuie\n}  — — ^ 


71-0 


Jc*0 


Applying  (D8)  and  the  equality 

n 

^2  n'  =  n-  (^+d 


Jt-0 


and  further  simplifying  (D12),  we  obtain 

Pr  {capture}  =  ^  Pr  {capture\n}  (n  +  l)  Gn(G+l) -(n+2) 

73  »  Q 


( D12 ) 


(D13) 


(D14) 
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